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Abstract. In this paper wc analyze American style of floating strike Asian call options belonging 
to the class of financial derivatives whose payoff diagram depends not only on the underlying 
asset price but also on the path average of underlying asset prices over some predetermined time 
interval. The mathematical model for the option price leads to a free boundary problem for a 
parabolic partial differential equation. Applying fixed domain transformation and transformation 
of variables we develop an efficient numerical algorithm based on a solution to a non-local parabolic 
partial differential equation for the transformed variable representing the synthesized portfolio. For 
various types of averaging methods we investigate the dependence of the early exercise boundary 
on model parameters. 

Key words. Option pricing, American-style Asian options, early exercise boundary, fixed domain 
transformation 

1. Introduction 

Asian path dependent options belong to the class of financial derivatives whose 
payoff diagram depends not only on the underlying asset price but also on the path 
average of underlying asset prices over some predetermined time interval. Such path 
dependent options can be often found at commodities markets such as oil, grain 
trade, etc. At expiration, the payoff diagram of such options is less sensitive with 
respect to sudden changes of the underlying asset value. Therefore a holder of an 
Asian option can effectively hedge the risk arising from a sudden price jump close to 
expiry. Typically, the payoff diagram of an Asian path dependent option depends on 
either arithmetic or geometric average of the spot price of the underlying asset. Such 
contingent claims can be used as a financial instrument for hedging highly volatile 
assets or goods. We refer the reader to references [IIlilllllEllllElEOlinilTOllllllIl 
discussing qualitative and quantitative aspects of pricing Asian path dependent 
options. 

In this paper, we focus on a special subclass of Asian options. Namely, we will 
investigate the so-called average strike Asian call options. At the time of expiry 
t = T, a holder of such an option contract has the right (but not obligation) to 
purchase the underlying asset for the strike price given as the path average of 
underlying asset prices. This means that the terminal payoff diagram for such an 
option has the form: V{S, A, T) = maxjiS* — A, 0}, where S = St '^s the spot price 
of the underlying asset, A = At is the path average of the asset prices St,t £ [0, T], 
over the time interval [0,T] and T > is the time of maturity. 

Concerning the method how the path averaged asset price A = At is calculated 
at a time i G [0, T], wc can distinguish the following methods of averaging of the 
path Su,u G [0, t]: 
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• arithmetic averaged options, where the average 



is given by 



(1) 




• weighted arithmetic averaged options, with the average AJ"" is given by 



and a is an exponential weight function a(^) = exp(— A^) with the averaging 
factor A > 0, 

• geometric averaged options, where the average A^ is given by 



In Figure [T] we plot two different sample paths of the underlying asset price (solid 
lines) and their arithmetic, geometric and weighted arithmetic path averages. In 
the case of weighted arithmetic averaging with a weight factor A > we can observe 
that the average A™" approaches the sample path St when A +oo. On the other 
hand, the weighted arithmetic average A^ approaches the arithmetic average when 



In this paper we are concerned with American style of Asian floating strike op- 
tions giving its holder the right to exercise it anytime before the obligatory expira- 
tion time t = T . Our main purpose is to provide a numerical quantitative analysis 
of the early exercise boundary position for a floating strike Asian call option by 
means of a solution to the transformed nonlocal parabolic partial differential equa- 
tion derived by Bokes and Sevcovic pQ [TB] for the case of arithmetic averaging. The 
main goal of this paper is to analyze the dependence and sensitivity of the early 
exercise boundary with respect to various model parameters. 

The paper is organized as follows. In the next section, we recall the partial 
differential equation for pricing floating strike Asian options. We also recall the 
method of dimension reduction of the equation. Next we discuss American style 
of Asian options and the early exercise boundary. A key tool for derivation of our 
numerical algorithm is transformation of the reduced equation to a flxed spatial 
domain. The resulting equation is a nonlocal parabolic PDE with an algebraic 
constraint between its solution and the free boundary position. In Section 3 we 
employ the method of finite difference approximation. We discrctize our nonlocal 
equation in space and time. To this end we make use of multiplicative operator 
splitting method. Results of numerical computations are presented in Section 4. 
We discuss the impact of the averaging method on the form of the early exercise 
boundary. We furthermore analyze the dependence of the early exercise boundary 
with respect to various model parameters. We furthermore make a comparison of 
its position for different averaging methods. 

2. Partial differential equation for pricing average strike Asian options 

In this section we recall the parabolic partial differential equation for pricing 
Asian options. Derivation of the pricing equation is based on standard assumptions 
made on stochastic behavior of the underlying asset price St,t G [0, T]. Henceforth, 
we will suppose that St follows a geometric Brownian motion, that is. 




(3) 




A ^ 0+. 



dSt = (r- - q)Stdt + aStdBt, 
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Figure 1 . Examples of evolution of underlying asset prices (solid lines) 
and their path averages (dashed lines). 

where r > is the risk-free interest rate, g > is the dividend yield on the 
asset, Bt,t G [0,T], is the standard Wiener process. Although the aforementioned 
assumption made on St has obvious deficiencies like constancy of the volatility a or 
normality of the distribution of asset returns, we adopt this assumption throughout 
the paper. 

Since the payoff diagram V{S, A, T) = maxjS — A, 0} depends on both the spot 
asset price S and the path average AaXt = T so does the price V of an Asian option 
for < t < T. It means y is a function depending not only on the underlying asset 
spot price S and time t but also on the average A of the underlying asset price over 
the interval [0,t], i.e. V — V{S,A^t). In order to derive the pricing equation for 
an option price V one has to calculate its differential dV over a time interval with 
an infinitesimal length dt. The price S = St well as its path average A = At are 
stochastic variables and so does the option price Vt = T^(S't, ^t, t). To calculate the 
differential dV we have to find a relationship between the differential of the path 
average At and the time t. Taking a differential of ^ we obtain the expression: 

dAt 1 St -At 



-s 

■Arithmetic 
■Geometric 
■Weigli.?i = 1 
■Weigli.X = 10 
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for the case of arithmetic average A'}. On the other hand, for the case of geometric 
averaging we obtain from 

* - 'hi Srdr + - In S't - * 



Af dt Jo ^ t ' t 

Similarly, for the exponentially weighted arithmetic average with the weight 
parameter A > we obtain 

dA'f" _ X{St - A]"") 



dt l-exp(-At)' 

In all averaging methods under consideration, we can conclude the following rela- 
tion between the differential dA as a function of the stochastic variable S and the 
differential dt. 



(4) 



dA^Afl 4,i 1 dt 



where the function / is given by 
X — 1 



(5) f{x,t) 



t ' 

A(a:- 1) 
f — cxp(— At) ' 

ln(a;) 



A' 

for arithmetic averaging, 
for exponentially weighted arithmetic averaging, 
for geometric averaging. 



It means that the differential dA is a stochastic variable with the leading order term 
of the order dt (see Kwok [T^], Dai [3], SevcovicfTB]). Applying Ito's lemma (cf. 
Kwok [12) for the function V = V{S, A, t) and taking into account @ we conclude 
the stochastic differential equation for the option price V in the form: 



dV 



dV 

'as 



dS' 



dV 
'dt 



2 dS^ 



dV 
'dA' 



An^,t 



dt. 



Under the assumption of perfect replicability of the market and nonexistence of 
arbitrage opportunities (cf. Kwok [12]). it can be shown that the risk neutral 
price V = V{S,A,t) of an Asian call option is a solution to the parabolic partial 
differential equation 



(6) 



s 



A' 



dV 



dA 



satisfying the payoff diagram corresponding to the averaged strike Asian option, 
i.e. 

V{S,A,T) ^ max{^- A,0}, 

where S,A>0, t€ (0,T). 

It is also well known that the PDE for average strike Asian options allows for 
a dimension reduction by introducing a new state variable x and the function W 
defined as follows: 



(7) 



S_ 



W{x,T)^jV{S,A,t), 
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where t = T — t. After straightforward computations we obtain a parabolic PDE 
for the function W{x, r): 

dW dW ^d'^W 

(8) ^ + [/(^' ^ " ^) " + '^l - + " ^ " ^ 

where re (0,r), x > 0. A solution W satisfies the initial condition: 

W{x,{)) =max{a;-l,0}. 

A solution W{x, r) to the aforementioned equation ([5]) is defined on the fixed spatial 
interval < a; < oo. The function V{S, A, t) given by V{S, A, t) = AW{S/A, T-t) 
corresponds to the price of an Asian floating strike path dependent option for the 
so-called European style of contracts for which the option expires exactly at the 
time t — T. 

2.1. American-style of Asian call options. In this paper we are concerned 
with American style of Asian options (cf. [Zl[3l[12]). In contrast to European style 
of options, American style options can be exercised at any time until the obligatory 
maturity time t = T. The holder of such an option has the right to exercise it 
or to keep it depending on the spot price of the underlying 5*4 at time t and its 
history {Su, Q < u <t] prior the time t. The boundary between "continuation" and 
"stopping" regions plays an important role in pricing American-style of options. It 
can be described by the function {Al^t) i— >■ SI = Sf{At,t), where is the so-called 
early exercise boundary (cf. [HlElElinilllllllllinilEl)- According to Kwok [H] , 
the set 

£ = {iS,A,t) G [0,oo) X [0,oo) X [0,T),V{S,A,t) = V{S,A,T)} 

is the exercise region. For the case of a call option, there exists a early exercise 
boundary function 5*/ = Sf{A,t) such that 

£ = {{S, A, t) g [0, oo) X [0, oo) X [0, T), S" > Sf{A, t)}. 

It means that = Sf{At,t), where At is the path average of underlying asset 
prices {S'u,0 < u < t}. We can search the early exercise boundary function in a 
separated form 

SfiA,t) = Axfit). 

For more details we refer the reader to Dai and Kwok [Sj or Bokes and Sevcovic 
[U [16]. The corresponding spatial domain for the reduced function W = W{x,t) 
satisfying ^ is therefore given by 

0<x<p{t), re (0,T), 

where p(r) = xf{T-T). From the continuity of ^(5, ^, i) at {Sf{A,t),A,t) we 
conclude that 

(9) ^{Sf{A,t),A,t)^l. 
It follows from the payoff diagram that 

(10) V{Sf{A,t),A,t) = Sf{A,t)-A, 

for A > and t e (0,T). In terms of the new state variable x, we conclude the 
following boundary conditions for the function W{x,t): 

dW 

(11) VF(0,r)=0, W{x,t)^x-1, — (x, t) = 1, at x = p(r), 

for T e (0, T). The initial condition for W{x,t) is 

(12) W"(a;,0) = max{x- 1,0}, Vx > 0. 
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Equation ([5]) and boundary conditions represent a free boundary problem, 
because the spatial domain < a; < p(r) depends on the unknown free boundary 
function p which is a part of a solution of the problem. 

2.2. Fixed domain transformation. Following ideas of transformation method- 
ology developed by Scvcovic in [TB] (see also [T] ) we introduce a new state variable 
^ and the transformed function 11 = n(^, r) defined as: 

/ p(t)\ dW 
(13) ^^InJ^^j, n(e,r)=W^(.T,r)-x — (x,r). 

After straightforward calculations we obtain that n(^, r) is a solution to the fol- 
lowing linear parabolic equation 



(14) 
where 



2 de 



6(e,r)n = 0, 



dr 



2 ' 



x—pe~ 



For details of derivation of p4)) we refer the reader to [TJ [16] . The initial condition 
for the solution n(^, 0) is: 



(15) 



n(^,o) 



0, 



C<lnp(0), 
C>lnp(0). 



The limiting value p(0) of the early exercise boundary at expiry r = (i.e. t = T) 
for the continuous arithmetic average type of an Asian option has been derived by 
Dai and Kwok [3]. For the geometric average it has been discovered by Wu in [18] 
(see also Detemple [H p. 69]). For arithmetic weighted averaged floating strike call 
option it has been derived by Bokes and Sevcovic in [1]. The value p(0) is given by: 



(16) p(0) = <^ 



p^O) 



p"'°(0) 



l + rT 



,1 



l + gT 
A + r(l-e-^^) 



pf(0) = max{i,l}, 



arithmetic averaging, 



weighted arith. averaging, 



geometric averaging. 



In the case of geometric averaging, the auxiliary number x entering the expression 
of /o(0) is a unique solution to the transcendental equation: 



(17) 



xqT — rT + ln(i) = 0. 



Next we recall the boundary conditions for a solution 11. Taking into account 
equations (|lip we end up with the Dirichlet boundary conditions: 



(18) 

Since ^ 



n(0,T) 



n(c»,T) = 0. 



and ^{p{t),t) = 1 we obtain ^{p{t),t) = at x = p{t). 
In the limit x — > p(r), assuming the continuity of n(^,T) up to the boundary 
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^ = 0, we obtain 



Passing to the limit x -> p(t) in equation (|5]) we obtain the fohowing algebraic 
constraint between p(t) and the solution n(^, r): 

-ir ~ q)p[r) - Y^(0' ^) + ^(^M " 1) + /(/'M' - ^) = 0- 

Therefore we obtain the following algebraic constraint equation between the free 
boundary position p(r) and the partial derivative 9^11(0, r): 

(19) (7p(T)-r-^ — (0,r) + /(p(r),r-r) =0. 

Notice, that this expression contains term ^(0, r) and therefore this is not suitable 
for numerical scheme, because the whole solution is sensitive of this term. Bokes 
and Sevcovic [1] suggested an equivalent form of (|T9| . They integrated equation 
(fT4| with respect to ^ <E (0,oo). Taking into account boundary conditions for 
n(C,r) and ^(oo,t) and using equality ([T^ they derived an ODE replacing 

the algebraic constraint between the free boundary position q{t) and ^(0,r) of 
the solution 11. It has the form: 



2 



(20) = ^{\np{T) + Ii{tT)ddj+qp{T)-q- 

+ (r-/(p(T)e-«,T-T))n(^,T)de 
Jo 

Remark 1. Denote by x" ^ {1 + rT)/{l + qT). Then x^qT - rT + ln(a;'') < 
x°'qT^rT+x'^ — 1 ~ provided thatr ^ q. Since the function x i— > xqT — rT -\-\t\{x) 
is increasing we have x"^ < x and consequently p°'{0) < p^(0). Similarly, for any 
\ > we have p'^"(0) < p'^{0)- In summary, we derived the following inequalities 
for the initial values of the early exercise boundary 

(21) p'^^{0)<p^{0)<p^O), if r^q. 
Clearly, p'^'^iO) = p^iO) = (0) = 1 m the case r = q. 

Remark 2. The solution W — W{x,t) can be easily calculated from the solution 
n and the free boundary position p. With regard to U3\) we obtain: 

^ {x-'W{x, r)) - -x-^n (In (p(r)/x) , r) . 

Taking into account the boundary condition W{x, r) ~ x — 1 at x ~ p{t) and 
integrating the above equation from x to p{t), we obtain 

H^(x,r) = --^ U(T)-l + y^ e«n(C,r)dCj- 



Then the option price can be calculated from equality ViJa^A^t) — AVK(S'/A, T — t). 
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3. Numerical algorithm 

In this section, we make use of tlie numerical algorithm proposed by Bokes 
and Sevcovic in [l] [16] for calculating of the early exercise boundary position p. 
We notice that the original algorithm was derived for arithmetic averaged floating 
strike Asian call option. In this paper, we generalize it for the case of geometric 
and exponentially weighted averaged floating strike options. 

The algorithm is based on a finite difference discretization in space and time 
variables. We restrict the spatial domain to a finite interval ^ € (0,-L), where 
1/ ;:g> 1 is sufficiently large. For practical purposes, it is sufficient to take L « 3. Let 
fc = ^ > is a time discretization step and /i = > is a spatial discretization 
step. We denote by II-' time discretization of n(^, Tj) and — p{Tj), where Tj — jk. 
By we denote full space-time approximation of the value Il{£,i,Tj). Then the 
Euler backward in time finite difference approximation of (|14p reads as follows: 



r) 



2 de 



where c' = c(tj) and c(t) = ^lnp(T) + r — g. We prescribe Dirichlet boundary 
conditions at ^ = and ^ = L for the function II-' . As for the initial condition, we 
use the vector 11° where 11° = n(<^i,0). Following [1], we make use of the operator 
splitting method to the above problem by introducing an auxiliary intermediate 
step n^^2 that splits the problem into two parts: 

• Convective part 



(22) 



W^2 -W- 



(23) 



Diffusive part 







df 

r + X— - f{x,T - t) 



A solution n-' 2 to equation can be approximated by the explicit solution to 
the transport equation 



0, 



for > ^i^d r G (rj_i,Tj] satisfying the initial condition n(^,rj_i) = ^(^) 
and the boundary condition 11(0, r) = — 1. After some computations (for further 
details see e.g. Bokes and Sevcovic [I]) we end up with the following solution: 



(24) 




if i^i = Ci ^ 

otherwise. 



ln£^-(r-g)fc>0, 



SENSITIVITY ANALYSIS OF THE EARLY EXERCISE BOUNDARY 



9 



In order to derive full space-time discretization scheme we make use of the central 
finite difference approximation of equation (j23|) . We obtain 



= 



2h 



Hence the vector W is a solution of a tridiagonal system of linear equations 

(25) a^nti+/3m^+ 7^1+1 =nr^ 
for i = 1, 2, . . . , n, where 

Boundary conditions for II-' are given by: IIq = —1, 11;^^ = 0, for j = 1,2, ... ,m. 
The initial condition for II" is given by equations ([T5|) and ([TE)) . In order to deter- 
mine the free boundary position we take equation (PH)) into account. Applying the 
forward finite difference approximation we obtain 

/>oo />oo 

(26) lnp^' = \np^-^+ n^'\£_)d^ - {^d^ 

Jo Jo 

f ^2 roo 

+fc U + — - q(P-^ - (r - fipi-'e-i,T- r,)) W i^)dC 

As for the approximation of the integral {^)d£^, we use the trapezoid quad- 

rature method. If we rewrite equations into the operator form then we obtain the 
following system of nonlinear algebraic equation for the unknown vector W and 
the free boundary position at the time Tj : 

(27) p^=J'(ff), W-i^Tip^), Aip^)Ui =W-K 

where T{p^) is a solution of the transport equation given by ((24)) . A{p') is a tridi- 
agonal matrix given by (P5|) and In ) is right side of equation (|26p . System of 
equations (P?]) can be effectively solved by means of successive iterations procedure. 
For j > 1, we set 11-'''' = n^~^ and p^'" = p-'^^. Then {p + l)-th approximation of 
II-' and p' is a solution of following system: 

(28) p^'^P+i = TiW'P), 

(29) n^-ip+i ^ W'P+^), 

(30) y^(pJ'P+i)nJ'P+i = 

We repeat the above iteration procedure for p = l,...,pmax until the prescribed 
tolerance \pj'P+'^ - p''P\ < toll is achieved (see Table[2]). 
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Input model parameters: g > 0, ? , <t, T > 0, A > 0, 



Input numerical parameters: n,m,L,pj^ 
Initialization: 

k = T/m 
h = L/n 



, toll 



end 



max{Y±^,l}, 
max {i, 1} , 

- A+r(l- 



for j = 1 to 



1, 



/ / arith. avg., 
// geom. avg., 
// weighted arith. avg., 



0, ^>lnp' 







Pi 

for p 



to Prnax- 

Solve: Aipj'P+^)W'P+^ = ff-iP+i for U^'P+^ 



if i\p 



endif 



,j,p+i 



< toll) 



break 



end 



Figure 2. A flowchart of the numerical algorithm. Input model 
parameters: r is the interest rate, q is the dividend yield, T expi- 
ration time, n is the number of spatial grid points, m is the time 
step, L = 3, A is the weight parameter. 



4. Numerical results 

The aim of this section is to present various computational examples of calcu- 
lation of the early exercise boundary position p for varying model parameters and 
averaging methods. In all examples discussed in this section we use the following 
numerical parameters: m = 10 000, n = 300, L = 3, Pmax = 500, toll = 10^*. In 
graphical plots we display the position of the early exercise p only. 

4.1. Free boundary for floating strike call options. 

4.1.1. Arithmetic averaged floating strike call options. In Figure[3]we com- 
pare the free boundary position p for various interest rates r = 0.02, 0.04, 0.06. 
Other model parameters are: T = 50, cr = 0-2, q = 0.04. We also compare the free 
boundary position computed by the algorithm described in Section 3 (blue line) 
and computational results obtained by Dai and Kwok in [3] (red dots). 

In Figure |4] we show the number of maximal inner iteration steps Pmax needed 
for achievement of the desired tolerance < toll ^ 1. We can observe that the 
algorithm requires considerably more inner iteration steps for small times < r <C 1 
for which the time derivative of pir) is large. As for the model parameters we chose: 
T 50, 0- = 0.2, q = 0.04, r = 0.06. 
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Figure 3. A comparison position of the free boundary position 
for various r — 0.02,0.04,0.06. We also compare our results with 
the method due to Dai and Kwok [3]. 
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Figure 4. A number of the inner- loop iteration in the algorithm 
needed to achieve prescribed tolerance. 



In Figure [5] we plot p(r) for different values of the volatility parameter a. Other 
parameters are: T = 50, r ~ 0.06, q = 0.04. It is worthwile noting that, for 0, 
the underlying stochastic process for St becomes deterministic. The option price 
at the time t can be therefore calculated as: St = S'oe'- Its arithmetic average 
is given by At = \Sq '^^ ^-q^^ • '"^^^ of g > r we have St < At for all 

t > 0. Therefore the corresponding option price is equal to zero because it is 
not worth to exercise it for any < t < T. In the case r = q, both St = At- 
Hence p{t) = 1. Finally, in the case r > g we obtain St > At, for all t > 0. 
Since we are dealing with American style of options we exercise the option at 
the time t* = arg maxtg[o,T] e-''(^-*)(S't - At). Denote by ^'(t) = ^^f^- Then 
t* = argmaxtg[o_T] 5'oe~'"-^e''*^'(t). Now it is easy to verify that 

^ ^ V r^^ir - qT^h- > 0, for ah t > 0. 
dt ^ in + 2 n! ^ ^' 

n— 
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^a=0.1 -0-a=0.2 -6-0=0.3-^0=0.4*0=0.5 




Figure 5. A comparison of the free boundary position p for large 
(top) and small (bottom) values of the volatility parameter a. 



Hence '^{t) is an increasing function. As both ^'(t) and e*"* we conclude that t 
In the case of an arithmetic averaged Asian call option we obtain from (H 
following explicit expression for the free boundary position: 



= T. 
the 



(31) 



p{t) = max < 1, 



1 + r(T - t) 

i + q{T- t) 



for a — 0. 



4.1.2. Geometric averaged floating strike caU option. In Figure|6]we com- 
pare the free boundary position p computed for the case of arithmetic and geometric 
averaging. Model parameters were chosen as: T = 50, cr = 0.2,(7 = 0.04, r = 0.06. 
Notice that p"(r) < p3{T) for ah r € [0,T]. 
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Figure 6. A comparison of the free boundary position for arith- 
metic and geometric averaging methods. 
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X 



Figure 7. A comparison of the free boundary position for expo- 
nentiaUy weighted arithmetic averaged Asian options for various 
weight parameters A > 0. 

4.1.3. Weighted arithmetic averaged floating strike call option. In Fig- 
ure [7] we plot the free boundary position p for various weight parameters A = 
0.001,0.1,0.2,0.5,1. 

It is easy to verify that lim = St- As a consequence we deduce that the 

A— >oo 

option price is equal to zero. This is also reason for the limit limA->oo Pa = 1- In 
what follows, we will estimate rate of convergence of pA — > 1 using the so-called 
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Table 1 . Experimental order of convergence a for the difference 
IIpa - PooIIoo = 0(A~") as A oo. 



A 


1 Pa — Poc \\oo 


Q 


0.2 


0.561828 




0.5 


0.413783 


0.333 


1.0 


0.320136 


0.370 


2.0 


0.247010 


0.374 


3.0 


0.212705 


0.368 


4.0 


0.191862 


0.358 


5.0 


0.177658 


0.344 


10.0 


0.147227 


0.271 
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experimental order of convergence. Assuming that 

\\p\ - Poolloo = ^max^lpA(r) - Poo(t)| = 0(A""), 

as A — > cxj, we can approximate the order parameter a as follows: 
(32) ln(A2) - In(Ai) ^ 

ln(||/3Ai - Poolloo) - ln(||/5A2 - Poolloo)' 

where Ai < A2. The results shown in Table [1] indicate a ~ 1/3. It means that it 
might be reasonable to conjecture that 

il/OA - Poolloo = 0(A"^/^), as A 00. 

4.2. Comparison of the early exercise boundary position for various av- 
eraging methods. In Remark 1 we pointed out that ^"""(0) < p°'{0) < p^{0) if 
r ^ g and /9"'"(0) = p"(0) = pS{0) = 1 in the case r = q, where p'',pS, p""" are the 
free boundary position for arithmetic, geometric and weighted arithmetic averaged 
Asian call options, respectively (sec (|2ip V This relation has been rigorously derived 
for T ~ only. However, it follows from results depicted in Figures |6] and [7] that 
it might be plausible to state the following conjecture on the comparison of early 
exercise boundaries for various averaging methods: 

Conjecture 1. Let p°',p^^p'"'^ be the free boundary positions for arithmetic, geo- 
metric and exponentially weighted arithmetic averaged Asian call options. Then for 
any < t <T it holds: 

p-^ir) < p^{t) < p^r). 

The rigorous proof of this conjecture based on the analysis of the governing 
equation (jl4p with constraint (|19l) remains an open problem. 



4.3. Initial early exercise boundary position. In this section, we investigate 
the initial early exercise boundary position Xf{0) at t = 0. It corresponds to the 
value g{T) = x/(0) at r = T. 

Let us denote by p(r) = p(r; r, q, cr, T) the free boundary position as a function 
of T G [0,r] and remaining model parameters, i.e. r,q,a^T. Using obvious scaling 
properties of the governing equation and the algebraic constraint ([T^ we can 
conclude that for arithmetic, geometric or weighted arithmetic average the following 
scaling property holds true: 

p{t; r, q, a,T) = p [rj^^Tr, Tq, VTct, 1^ 
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Figure 8. Dependence of p{T) on the parameter r (top), q (mid- 
dle) and tr (bottom). 



With regard to the previous argument, it is therefore sufficient to study dependence 
of p{T) on r, q, a for arbitrary but fixed value of the parameter T . In Figure[5] (top) 
we plot dependence of p{T;r,q,a,T) on r e (0.005,0.1) for a = 0.2, q = 0.04. In 
Figure [5] (middle) we plot dependence of p{T;r,q,a,T) on q & (0.05,0.1) for a = 
0.2, r = 0.06. Finally, in Figure |8] (bottom) we plot dependence of p{T;r,q,a,T) 
on a e (0.01, 0.4) for r = 0.06, q = 0.04. In all case we chose T = 50. 

We furthermore proposed and consequently estimated the analytical formula for 
dependence of p(T; r, a, T) on r, tr, T. With respect to the scaling property we 
can fix T = 50. Among various nonlinear estimators for the value /5(T; r, cr, T) on 
r, g, (T, 1 the best estimation results we achieved with the function 

e''PP(r; r, g, a, T) ^ g'^^^iX- Tr, Tq, VTa, 1) 1 + ( + -P^ 

\Pir + P2qJ q 

with /?i = -0.15064, = 7.74793, = 0.79067, /34 = 0.09193 with RSS ^ 
8.6723.10^4 where as an error indicator we chose the least square minimizer: 

M 

RSS ^ Y.^g'^'"' ' * ' "1 ' ^) - (^))' ' 
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where T ^50 and parameter samples rj,qj,aj e [0.01, 0.11] x [0.01, 0.11] x [0.2, 0.8] 
were generated from M = 100 random vectors. We denoted by Pj{T) the nu- 
merically computed free boundary position p at r = T computed for the model 
parameters r = rj , q ~ qj , a — Uj and time horizon T . 
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